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1 Introduction 



SU{N) group integration formulas at large N have attracted interest since 
the late 70's. They occur in the study of large N lattice gauge theories and 
in connection with matrix models. The latter have found a wide range of 
applications over the years in areas going from 2-dimensional quantum grav- 
ity to 4-dimensional Yang- Mills theory. The class of matrix models involving 
large group integrations are the unitary matrix models (see Ref. |l|] for a 
recent review). 

Within lattice gauge theories, SU{N) group integrations arise naturally 
when studying strong coupling expansions of lattice gauge theory models, 
which simplify in the large limit. The most standard method uses an 
effective action approach centered around the one-link integral. Our motiva- 
tion for this paper arose in the context of A/" = 1 SUSY Yang-Mills on the 
lattice, an interesting testing ground for non-perturbative methods, and for 
which some Monte Carlo data are available [0]. The model contains lattice 
gauge fields and quarks (gluinos) in the adjoint representation of the group. 
When studying the model at strong coupling one notices that the effective 
potential method cannot be carried through, lacking an equivalent one-link 
integral formula in the adjoint representation. Our approach Q was based 
on the hopping parameter expansion expressed as a sum of random walks 
on the lattice |^ ^f. Then, one has to perform the gauge field integration 
of an expression containing products of traces of Wilson loops in the adjoint 
representation. It turns out, as will be shown in this paper, that the leading 
large contribution is rather trivial from the point of view of large in- 
tegration. It arises from loops which are pure-spikes (or pure backtrackers) 
and that actually do not depend on the gauge fields. 
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Another instance in which group integrals of quantities in the adjoint 
representation have acquired relevance is in the context of induced QCD 
The Kazakov-Migdal model has scalar quarks in the adjoint representation 
coupled to a /5 = lattice gauge field. A wide variety of techniques have 
been applied to this model. The strong coupling (large mass) region solution 
of the gaussian model was obtained by Gross Again, as realised later by 
Khoklachev and Makeenko 0, the role of the large N integration over the 
gauge variables is simply to restrict the sum over random walks to pure-spike 
paths. Their contribution does not depend on the gauge fields, so that it can 
be taken outside the integral. 

In this paper we will address the problem of computing the leading non- 
trivial large N integral of a product of traces of Wilson loops in the adjoint 
representation. The result is more easily presented in the form of the evalu- 
ation of a generating function (Eq. (|1])) for adjoint Wilson loops. The main 
result is presented in Eq. (|T3|). The formula is obtained by expanding the 
exponential, performing the integration and resumming the result, so that it 
is only valid for sufficiently small lattice couplings 5(7), i.e. it is the strong 
coupling expression of the generating functional. The result is the exponen- 
tial of a quantity which is of order 1 (~ N^). This contrasts with the typical 
A^^ behaviour for the fundamental representation. Indeed, in applications to 
strong coupling expansions of theories with matter fields, there is an addi- 
tional A^^ — 1 term arising from trivial Wilson loops which are proportional 
to the trace of the unit matrix in this representation. This term comes from 
the pure-spike paths mentioned previously. Thus, the terms we are calcu- 
lating here give large A^ subleading contributions to the observables of these 
theories. We emphasise nevertheless that our calculation gives the leading 
term in which the actual details of the group integration play any role. 
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The paper is organised as follows. The main quantity we are calculating 
is given in Eq. (|I|) of Section 2, and the result for it is given in Eq. (0). The 
rest of the section contains the outline of the derivation, with details and 
terminology collected in the two appendices (this includes a useful diagram- 
matic formalism to deal with the computation of group invariants). Readers 
not interested in the derivation might jump directly to Section 3, where the 
formula in question is put to work in calculating large corrections to the 
strong coupling values of observables in SU{N)/Zn lattice gauge theories, 
with or without matter fields. Explicit formulas are given for the lattice ver- 
sion of A/" = 1 Supersymmetric Yang-Mills theory. This includes the large 

corrections to the gluino condensate and the leading non-OZI contribu- 
tions to the scalar meson propagator. Similar calculations can be carried on 
for the gaussian Kazakov-Migdal model. As an example, we show how from 
our formulas one can deduce the leading behaviour of the propagators of 
the p-scalar bound states. The masses of these states are explicitly given. 
Finally, the paper is ended with a short conclusion section. 



2 Group integration formulas 

We consider an SU{N) gauge theory defined on a hypercubic, d-dimensional 
infinite lattice A. To any link i of the lattice there is associated a group 
element ge. Most expressions that we will use involve the matrices corre- 
sponding to the adjoint representation of the group UA{ge) = Ug. Regard- 
ing geometrical objects on the lattice, such as lattice paths, we will adopt 
the terminology and conventions of Ref. |@]. Given a closed lattice path 7, 
Ti{Wa{i)) labels the trace of the corresponding Wilson loop evaluated in 
the adjoint representation. 
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Our main goal in this section is to compute, to leading order in N, the 
result of integrating over the group variables an arbitrary product of traces 
of adjoint Wilson loops. This result is best collected in the form of the 
following generating function: 



In this expression, d/ijy(g'^) is the Haar measure on the group variable g^, and 
the product extends to all the links £ of the lattice. The sum in the exponent 
runs over all equivalence classes of closed paths. Each equivalence class 7 
contains all paths such that the traces of the adjoint Wilson loops computed 
along them coincide for an arbitrary gauge field configuration. The common 
value of the trace of the adjoint Wilson loop is labelled Tr(VFyi(7)). The 
arguments of the generating function are the real parameters s{j). Notice 
that Z{s) can also be looked at as the partition function of a certain kind 
of pure gauge theory. The parameters 5(7) would in that case appear as 
the (lattice) couphngs of the corresponding action. Although the action is 
not the most general SU{N)/Zn- invariant action, other terms would give 
superfluous information. From this point of view the result of this section 
can be interpreted as the evaluation of the partition function of this pure 
gauge model in the strong coupling region (small values of 5(7)). 

The expansion of the exponential in Z gives us the integrals of traces that 
we are interested in: 



Eili...%) = ll{ dii4ge))T:t{WAiji))...Tr{WAijn)) ■ (3) 




(1) 





with 
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To compute the previous expression one must first select a representative of 
the class 7. We will choose a closed path without spikes (backtracking parts). 
It is always possible to reduce an arbitrary path by eliminating its spikes in a 
systematic way and reach a path which is called simple in Ref. [H. As we will 
see later there is more than one simple path within each class, but anyone 
will do. Let us, from now on, use the same symbol 7 to refer to the class and 
to its representative simple path. In the rest of this section we will show the 
main steps in the evaluation of expression @. Many of the technical details 
are collected in Appendices A and B. 

A possible approach would be to express the adjoint Wilson loops in 
terms of fundamental representation ones and then use the formulas available 
for the latter. Instead we would use a procedure which could allow the 
calculation of traces of adjoint Wilson loops and other gauge invariant objects 
systematically to all orders in What we do first is to express the 

integrand in Eq. (Q) as a sum of products of the elements of the adjoint 
representation matrices Ue: 

The sum extends over the 2L{N'^ — 1) possible values of the indices. For 
simplicity we will refer to the indices jointly by the set X = {oi, 61, . . . ai, 
The integer L gives the length of the total joint path 71 . . . 7„. The labels 
ii . . .i^ refer to the links of this joint path taken in an arbitrary order and 
where a given link can occur more than once in the sequence. Finally, the 
symbol 6{p) contains the contraction of the indices necessary to produce the 
traces given in Eq. (^). It is given by a product of L Kronecker deltas. It 
is clear that which index is contracted with which depends on the naming 
of the indices and on the geometry of the paths. This information can be 
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coded into an application p : T — T which maps the contracted indices onto 
each other (p(c) = d ; p{c') = c). Since the apphcation is bijective it can be 
seen as an element of the group S2L of permutations of 2L indices, having 
the properties of a pairing. It is clear from the definition of the index set X, 
that there is another pairing r G S2L naturally defined on our system. It is 
the application which maps the two indices of a given link matrix onto each 
other: r(aj) = 6j ; T(6j) = Oj. 

The next step in the calculation is the evaluation of the group integrals 
involved in Eq. (^). Our procedure will be to express the adjoint representa- 
tion matrices in terms of the fundamental representation ones, and then to 
make use of the available formulas for that case. Here we will just explain 
and make use of the results, collecting the details of the derivation within 
Appendix A. First, one can look at the result of integrating over a single link 
variable. The result can be expressed as a linear combination of invariant 
group tensors, with coefficients depending on A^. These tensors have indices 
within those elements of I which belong to the matrices corresponding to 
this link. The result of integrating over all link variables can be expressed 
similarly as a linear combination of invariant tensors, this time involving all 
the indices of the set X. Thus we can write: 

£(71 ...in) = J2 ^{pr E ^^(^) n^r , (5) 

where X('(9)'^1"'^2l jg invariant tensor and Ci\f{'d) the corresponding A^- 
dependent coefficient. The indices q simply label the elements of X. The 
symbol stands for a permutation {i) : I I) . In Appendix A (Eq. we 
describe how one can associate to any permutation of a certain type an in- 
variant tensor ("i? — > T{d)). There, one can also see that if the permutation 
has a nonvanishing number of 1-cycles, the corresponding tensor vanishes. 
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It is clear by the construction that both the invariant tensor and its coeffi- 
cient factorise into a product of contributions of single links. This and other 
properties (like the absence of 1-cycles) restrict the space of possible The 
full set labelled Sr defines a certain subgroup of 821- We will leave out its 
characterisation for later and continue with the outline of the procedure. 

The final step consists in performing the summation over the values of the 
indices, which is equivalent to contracting the colour indices of the invariant 
tensor with the S{p). Defining 

D^i^, p) = J2 Sipr T(^9)=--- (6) 

we finally get: 

^(7i---7n)= 5^C;v(^)/^^(t9,p) . (7) 

The result is expressed in terms of the coefficients Ciy{'d) and Diy{{},p). 
Explicit forms for these quantities can be traced back from Eqs. ( PT| ) and 
(]4^ ) in Appendix A. The main point now is to analyse the behaviour of these 
coefficients as ^ 00. 

The large- A^ behaviour of the coefficients Cjv('i9) follows from the result 
for the corresponding coefficient in the one-link integral, Eq. (P^D; one gets: 



Cn{i^) = ^(^9) y—j + subleading terms . (8) 

where the leading prefactor A{'d) is zero for all but a restricted class of 
permutations ^9. It is not hard to characterise this subset within Sr- It is 
given by those permutations satisfying: 'dr'&T = e, where e is the identity 
mapping. For these permutations the coefficient A{'d) is equal to 1. This 



result rests on the work of Ref. |T0 |. 

The machinery to deal with the contractions of invariant tensors entering 
the definition of Dn, and extract their leading A^ behaviour, is set up in 
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Appendix B, whither we refer for details. The method follows by setting up 
a Feynman- diagram interpretation of each term in the sum Eq. (|^. For that 
we need to enter into the structure of the invariant tensors Tid). These are 
products of traces of products of the group generators in the fundamental 
representation (cf. Eq. To any p-cycle contained in the permutation d 

there corresponds an invariant tensor of the form Tr (A'^^ . . . A^'') (where A^* is 
a group generator). It is clear why the tensor vanishes if there are 1-cycles. 
We look at these traces as vertices of a Feynman diagram where the indices 
play the role of the legs. The vertices are not symmetric with respect to 
the exchange of the legs but only cyclic-symmetric. This form of vertices 
is characteristic of planar diagrams in large hermitian matrix models, for 
example. Thus, for the full d one obtains a bunch of vertices with a total of 
2L legs. The contractions encoded in p can be seen as propagators joining 
pairs of legs of different or the same vertices. Thus to every term in Eq. (|^ 
one can associate a (vacuum) Feynman diagram with V vertices and L lines 
(note that L, which is the total number of pairs of colour indices in Eq. (^, 
is the only topological quantity of the diagrams being equal for all the d and 
p). We further call K the number of connected components of the diagram. 

In Appendix B we prove that the large N maximal leading contribution of 
a diagram with K connected components, V vertices, and L lines is given by 
m2~^N'^^~'^~^^ , where m is an integer that can be zero or one. Corrections 
go like inverse powers of A^^. If we combine this with (|]), we obtain that the 
leading term coming from CN{^)DN{i), p) goes at most as N"^^'^ . 

Now, from the topological constraints it is clear that 1 < K <V . There- 
fore, the maximum value of 2K — V would occur when K = V . This means 
that all propagators go from one vertex to itself. To understand the impli- 
cations of this we need to make use of some of the properties of the vertices 
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appearing in such a diagram and its relation to the original set of lattice 
paths. Since vertices arise from link integration they can only contain in- 
dices pertaining to matrices of the same link. Furthermore, vertices only mix 
the first indices of the matrices among themselves and the second among 
themselves (see Appendix A for details). Hence, a contraction between two 
indices of the same vertex comes from a contraction of the link matrix with 
itself. But since the first index is contracted with the first or the second with 
the second, we have either UllAi or UdAl, which due to the orthogonal nature 
of the adjoint representation matrix allows to eliminate the two matrices al- 
together. In terms of the original paths a self contraction corresponds to 
the path moving in one direction forward and backward consecutively. This 
gives what we call spikes. But remember that our original paths were sim- 
ple, and hence had no spikes. In conclusion, the restriction to simple paths 
tells us that there are no self-contractions of our diagrammatic vertices, cor- 
responding to the situation when these vertices come from normal ordered 
operators. 

We then conclude that all connected components of the diagram have at 
least two vertices. Hence, the leading term arises from diagrams for which 
2K = V and which give a contribution of order one as — > oo. Finally, we 
arrive at 



where uo is an integer summing the contribution from all the leading order 
diagrams. This is the main result that was needed in Ref. to justify the 
restriction to pure-spike paths. In what follows we will proceed to compute 
the value of z/q. 

First of all let us study what can we learn about the graphs giving nonva- 
nishing leading contribution {0{N^)). We have not yet given the condition 




(9) 
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that has to be satisfied by ■& to give a nonvanishing leading contribution to 
Dn{'&,p)- This occurs (see Appendix B) when the permutation ■& satisfies 
the relation: dpiSp = e. This is very similar to the characterisation coming 
from the hnk integration, with the replacement of r by p. Our final number 
z/q comes from counting the number of possible permutations d within Sy 
satisfying both relations simultaneously. Furthermore, the set Sy is charac- 
terised by saying that it contains all permutations having no 1-cycles and 
such that all cycles contain indices belonging to matrices of the same link 
and to either the first or the second index of these matrices (but not both) . 

To understand the implications of the previous conditions, let us return 
to the connection between the original paths 7j and the permutation p. Any 
of the paths is an ordered sequence of matrices. Since the matrices have two 
indices we get a classification of indices into two sets: those appearing as a 
first or a second index of the matrices. Notice that if the link is followed in 
the reverse direction the definition of what index is first or last is exchanged. 
Any path is then associated to two equal length closed sequences of indices. 
It corresponds to two of the cycles of the permutation pr. The two cycles can 
be connected by means of r (or p) . Now, from the conditions characterising 
the leading we conclude: 

pr = ^pr^-^ . (10) 

This means that •& must map complete cycles of pr into equal length cycles. 
Given the identification of these cycles with the original paths, we see that 
a leading '& just maps complete paths of equal length among themselves. 
But since '& only connects indices belonging to the same lattice point, we 
conclude that the actual paths connected by ■& have to coincide. There are 
a few very simple conclusions one can extract from this result. The first is 
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that there is a factorisation property satisfied by the integral £^(71, . . . 7„) . Its 
value factorises into a product of contributions of distint paths. This follows 
from a similar factorisation property for 'd, since indices from different paths 
cannot be connected by it, despite the fact that the paths might overlap 
partially. Thus, the final number of maximal permutations Uq factorises as 
well. One might have guessed this factorisation from the standard arguments 
in the large N limit; notice, however, that the factorisation is not complete: 
as we will see, the integral of a given Wilson loop to the power n is not the 
n-th power of the integral of a single Wilson loop, so that the result is not 
classical. 

It only remains, therefore, to compute uq for the case of a single Wilson 
loop Tr {Wa{^)) taken to the power n. To do this, it is convenient to express 
the path 7 in the form (7')^^, where c<j is a positive integer. This notation 
means that 7 is built up by describing 7' sequentially lu times; 7' itself must 
not be reducible equivalently. We will refer to cu as the winding number of 
the path, and write it as a; (7). These considerations are necessary, as our 
results will depend explicitly on a; (7). 

For our path configuration we have that pr has n cycles of length ^(7) 
(the length of the path). The permutation maps one cycle into other 
preserving the ordering within each cycle. Thus, it is specified by giving the 
corresponding mapping of the indices at a given space point. Naively one 
would say that the number is then n\. In fact this has to be corrected by two 
facts: the first is that we have to eliminate all permutations -d having 1-cycles; 
the second is that one must take into account that if a; (7) is larger than 1, 
the mapping has to respect the cyclic ordering in which 7 describes these 
points. This last issue can be dealt with by arranging the indices belonging 
to the same path into cycles of length a; (7). In a compact mathematical 
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form we have that, given a permutation t' having n cu-cycles, we want to 
count how many permutations d' without 1 cycles of the nuj indices satisfy: 
d'r'd'^^ = t'. Let us call this number PJ. The action of is simply the 
mathematical statement of a change of names of the indices, so that without 
the requirement of no 1-cycles, the solution would be simply n!u;". This 
comes because we must map complete cycles into complete cycles, so we 
have n\ possibilities. Then after having selected the mapping of the cycles, 
there are u possible assignements for each cycle (it is enough to specify the 
mapping of 1 element of the cycle and the rest would follow): this gives the 
w". Subtracting the permutations which have 1-cycles takes a little more 
effort but presents no real problem. As a matter of fact, it is easier and 
more compact to determine the expression for the generating function first 
and deduce from it the expression of the quantity P". For the generating 
function corresponding to a single path one obtains: 



For completeness we give an expression for the numbers PJ^ themselves: 



Finally, we combine the result Eq. (|T^) for all the paths via the factorisa- 
tion property and arrive at the desired leading contribution to the generating 
functional Eq. i^: 



which is the main result of this section. As explained in the introduction the 
exponent is of order A^*^. Had we added a contribution to the initial action 




(11) 




(12) 




(13) 
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proportional to the trace of the unit matrix in the adjoint representation, the 
term would have remained unchanged by the integration (being a constant) 
and would have given a contribution to the exponent of order A^^ — 1 (the 
dimension of the adjoint representation). Although trivial, this remark is im- 
portant in practice because such a term actually arises in many applications. 
It can be associated to the trivial path of zero length. 

3 Application to strong coupling expansions 

In this section we will use the group integration formula obtained in the last 
section to derive results which are relevant for the behaviour of SU{N) /Z^ 
gauge theories on the lattice, with or without matter fields, at strong cou- 
pling and large A^. These include pure gauge theories, gauge theories with 
fermions in the adjoint representation, a particularly interesting case being 
the lattice version of A/" = 1 Supersymmetric Yang-Mills 0, and also gauge- 
Higgs models with scalars in the adjoint representation like the Kazakov- 
Migdal model 0. In all cases, after integrating out the matter fields via a 
hopping parameter expansion, one ends up with an effective action of the 
type given in Eq. (|1]). Our integration formula Eq. (p!3| ) allows us to write 
immediately the expression for the free energy in the strong coupling region 
at large A^. In what follows we will give the explicit formulas in some of the 
cases that have been studied previously in the literature. 

Consider first the case of SU{N)/Z]\r pure gauge lattice theories. For this 
case the consequence of our formulas can be easily deduced. As mentioned 
at the beginning of Section 2, the generating function Z{s) can be seen as 
the partition function of a pure gauge model where the couplings /3a(7) are 
identified with 5(7). The free energy of the model is then given by the 
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exponent of the right hand side of Eq. (T^). As a particular example we can 
write down the case of the Wilson-type plaquette action model. The free 
energy per unit volume of the model is given by: 

F{(3A)/Volume='^^^^{-(3A-\og{l- (3a)) , (14) 

which has a singularity at = 1. Actually, the position of the singularity 
coincides with the position of the third-order phase transition point in the 
2-dimensional lattice gauge theory [jll|. In higher dimensions the model is 



known to have a first order phase transition line at /?a < 1 [Hll, [TT[. Thus, 
the strong coupling phase depicted by Eq. (0) extends into the metastable 
region as is usually the case for these transitions. There is one atypical 
behaviour of the free energy Eq. ( p!4D which explains why this formula was 
not found before. This is the fact that it is of order 1 = 0{N^). Typically 
the free energy goes like iV^, as for example is the case in this same model 
above the first order phase transition point. The different N dependence of 
the free energy above and below the transition suggests that the position of 
the transition is given by the vanishing of the weak coupling free energy. This 
is indeed the case in 2 dimensions. An A^^ behaviour is also shown by the 
free energy of the mixed fundamental-adjoint action for any non-zero value 



of the fundamental coupling . Thus, in older work the free energy (|T^) of 
this model is taken to be zero Remarkably, from Eq. (0) one can easily 
read out the free energy for an arbitrary SU{N)/Zn invariant action. 

Let us now consider the addition of matter fields in the adjoint represen- 
tation. The general features and behaviour is similar for scalar or fermionic 
matter. For definiteness we will illustrate the formulas restricting to the 
fermionic case. Indeed, our concern with A/" = 1 SUSY Yang-Mills triggered 
our interest in the problem. Here we will follow the procedure of Ref. H. 
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Performing the integration over the quark (gluino) fields leads to a deter- 
minant. In the large mass (small hopping) region, this determinant can be 
written as the exponential of an effective action expressable as a sum over 
closed paths on the lattice: 

oo ^ 

^e// = ^/EE E yTr(W(7))Tr(r(7)) , (15) 

xeA 1=1 -y€Siix) 

where Si{x) is the set of closed paths of length / with origin in the lattice 
point X, and r(7) is a spin matrix associated to the path, whose explicit 
form we will detail later (cf. Eq. ([I8|)). The factor rjf takes into account 
the different possibilities for the nature of the fermion field: rjMajorana = 1/2, 

while rjDirac = 1- 

In Ref. it is shown how it is convenient and possible to rearrange the 
summation over paths into a sum over classes of paths labelled by what is 
called a simple path. A simple path is a closed contour that has no spikes or 
backtracking pieces. It is easy to see that for all paths associated with the 
same simple path the trace of the Wilson loop is exactly the same, so that 
they actually belong to a single equivalence class of the type defined at the 
beginning of Section 2. However, the classes considered in this paper contain 
more than one simple path and its corresponding classmates. The reason is 
that in Ref. closed paths are defined as ordered sets of links, which have 
therefore a first and a last element. It is clear, however, that which link is 
taken as the first one is irrelevant for the gauge field dependence, being given 
by a trace. This is just a technical point, but essential to get the numbers 
right. To compensate for this point one has simply to multiply by a counting 
factor which gives the number of links that can be chosen as first element. 
Curiously this factor is not the length of the path but l(^)/u(^), where 
uj{^) is the winding number factor defined in the last section. When the path 
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goes round a given loop uj{^) times, it is clear that the relevant factor is the 
number of links in the basic loop, and not in the total path. In addition, the 
orientation of the path is irrelevant for the gauge field dependence, so that 
our present classes contain each path and its reverse. In summary, we stress 
that the classes of paths defined in Section 2 contain a basic simple path, 
its reverse and all the /(7)/u;(7) simple paths obtained by selecting a point 
along the path (as origin). 

Finally, one can rewrite the effective action (Eq. ([T5|) ) in the form of 
the exponent of Eq. (|T|), with coefficients 5(7) which depend on the basic 
parameters of the original theory. These coefficients are 0{1) in N. One can 
include in them the contribution of the pure gauge action. Restricting for 
definiteness to one flavour adjoint fermions at zero gauge couplings, which 
contains the lattice version of A/" = 1 supersymmetric Yang-Mills, one obtains 
the following coefficients: 

^(7) = ^ , (16) 
r{l) = (Tr(r(7)) + Tr(r(7-^)) , (17) 

where 7"^ is the reverse of the path 7, and the matrix r(7) is the spin matrix 
of the path, given by the product along the path of the link matrices: 

/€(rI,-7„) . (18) 

In the previous expression k is the hopping parameter, r the Wilson parame- 
ter, and the spin matrices 7q, are just the Dirac matrices up to a sign, positive 
if the path traverses the link with positive orientation and negative if not. 
The factor 1 — ^ in Eq. ([T7|) is the result of resumming all paths within each 
class, with ^ given by: 



1 - VI -4(2rf- l)/t2(r2 - 1) ^ ^ 

e = — ^ — 2 ' ^^^^ 
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where d is the space-time dimension. 

With the previous expression and formula ([I3|), one can obtain the free 
energy of the model at large N: 

F = r]f {N^ - 1) Fq Volume - 

where the sum runs over the equivalence classes of paths 7 defined in Sec- 
tion 2. 

We see that our result provides the 0{1) contribution to the free energy. 
This is not however the leading piece, which is given by the term proportional 
to Fq, and is of order N"^. The latter was given in Ref. |^ and its value is: 

Fo = r/^(-dlog(l-0 + (rf-l)log(l-^|^0)Tr(I.) , (21) 

where is the identity matrix in spin space. This term arises from the 
paths whose gauge field contribution is proportional to the trace of the unit 
matrix. This occurs for the so-called pure-spike paths, which are paths such 
that after eliminating all backtracking parts we are left with a single point 
(a path of zero length). This contribution is, however, fairly trivial from the 
point of view of the group integration, so that the 0{1) contribution is the 
first one in which the characteristics of group integration show up. 

Obtaining a closed expression for the 0{1) contribution to the free energy 
involves performing the sum over path classes explicitly. Similar resumma- 
tions can be performed (see Ref. for a derivation of these resummations 
and for earlier references on the subject). The appearance of the factor ^7(7) 
in the formula has prevented us so far from obtaining such a closed expres- 
sion. 
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Other interesting quantities can be computed to the order in at which 
we are working. This includes the gluino condensate, which can be obtained 
by an appropriate differentiation of the free energy. As in the previous case, 
the leading term given in Ref. [Q comes from pure-spike paths, and is (up 
to some multiplication factors) independent of the group representation in 
which the quarks are living. This is a particular manifestation of the triviality 
from the point of view of group integration of the leading term. Finally, our 
expression is as follows: 

(Ar2-l)Tr(I,) 



(^(x)^(x)) 



1 y-Ml^(^^ + ^^] (22) 



Volume ^ rjf (1 — 2,^) \ 1^(7) 1 — ^(7) 



7 

where ^(7) is the length of any representative simple path. The division 
by the volume in the previous expression is simply accounted for by picking 
a single simple path 7 from the space of all space-time translates. The 
correction is the first unquenched contribution. The leading one coincides 
with the result of the quenched approximation [Q. To arrive at expression 
Eq. ( p2D we have made use of the resummation formulas of Ref. ||^. The 
origin of the 1 — 2,^ factor can be traced back to the rearrangement of sum 
over paths into sums over simple paths. 

In a similar fashion, it is possible to derive the large N correction to the 
meson and p-gluino propagators and spectrum, which are given in Ref. 
We will not give the formulas explicitly since they are notationally fairly 
lengthy. Instead we will provide the expression of a quantity for which our 
present computation gives the leading large behaviour. This quantity 
is the non-OZI contribution to the meson propagator. In a theory with 
different equal mass flavours, this term gives the main contribution to the 
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mass difference between the non-singlet and singlet mesons. It contains what 
in the continuum would be the contribution of the anomaly. When expanding 
the fermion propagators in terms of paths, the non-OZI contributions come 
from 2 closed paths originating at the 2 meson positions (see Ref. 0] for 
details). Due to the main factorisation property which we have obtained for 
Wilson loop integration at this order in A^, the contribution is only non-zero 
when the two paths (classes) coincide. To the order in at which we are 
working, this term is non-zero both in the quenched approximation and in the 
full theory. For simplicity, we give the contribution to the scalar propagator: 




where the sum now extends over all closed terse paths 71 with origin in x 
and 72 with origin in y (see Ref. Q for an explanation of the lattice path 
terminology and an account of the resummation formulas). The S^-^^y.^ imposes 
that the two paths have the same common reduced simple path, which we 
label 7. It is possible to eliminate the constraint given by the 6 in the previous 
formula, but its treatment exceeds the purpose of this paper. 

In a similar way one can study other models, like the Kazakov-Migdal 
model which has scalars in the adjoint representation. The gaussian model 
can be treated in a similar way as the fermionic model. In this model there 
are composites formed by p scalars, which are the analog of the p-gluino 
bound states of the fermionic model The operators which produce these 
states can be written in our notation as follows: 

0,{x)^<PaA^)---<P.^{x)T{l^Y-''' ■ (24) 

At strong coupling the correlator of two operators at different points (say x 
and y) can be written as a sum over p paths (71 . . . 7p) joining the fields at 
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the two different points. The gauge field dependence of such a contribution 
is as follows (we restrict to correlators that are leading at large A^): 

^(|_(p)^)ai...ap ^(|_(p))fep...fei ^Wa{1i)T'''' ■ ■ ■ (WAilp))""'''' • (25) 

Now we have to evaluate its expectation value with respect to the gauge 
field in the large N limit. The quantity is not a product of traces of Wilson 
loops. To handle this in a simple way with our formulas, we might add a 
new big link I going from y to x and its corresponding group element (as an 
independent variable). Now, one can form p closed paths by attaching to each 
of the original open paths 7^ the big link / (note they all have u = 1). One 
can then compute the expectation value of the product of the corresponding 
Wilson loops. Clearly s for these paths is zero, since none of the paths 
in the effective action contains the big link. The large leading term of 
the connected expectation value is obtained for all paths belonging to the 
same class. The corresponding contribution is the p-th derivative of the free 
energy in Eq. (0) at s = 0, which is {p — What has this to do with 
the expectation value of Eq. (pS])? To see that, one has simply to integrate 
over the big link group variable. One gets terms like Eq. (p5|) multiplied by 
{jj'Y- Then one has to multiply by the number of different p-cycles, which 
is (p — 1)!. Hence, finally, the expectation value of Eq. ( pSf ) to leading order 
is given by (y)^, and is achieved when all the path classes 71 . . .7p coincide. 
Similar results can be obtained in the fermionic model, thus justifying the 
treatment of p-gluinos of Ref. [Q. As an application of this result one can 
proceed in a similar fashion as in Ref. |^ and deduce the values of the masses 
Mp of the p-field bound states for this theory: 

eo.,M,^(ai^((l^)^^(^)-)-.,,., 
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a being an integer number in { — (rf — 1), . . . , (c? — 1)}, which accounts for all 
the possible doubler modes (see [^] for details), and ^ is given by the same 
formula as before (Eq. |19D with — 1 set to 1. 

4 Conclusions 

In this paper we have studied the leading behaviour for large number of 
colours of SU{N) group link integrals of products of traces of Wilson loops 
in the adjoint representation. We have proven that the leading (non-trivial) 
result is 0{1) = 0{N^). The corresponding coefficient has been computed in 
Section 2. The result factorises into contributions of Wilson loops associated 
to inequivalent paths. This allows us to compute up to this order and in the 
strong coupling (small 3(7)) region the corresponding generating function, 
given in Eq. (|T3|) of Section 2. 

In Section 3 we have shown how this generating functional and its deriva- 
tives enter in expressions which are relevant in the context of strong coupling 
expansions for gauge theories with matter fields in the adjoint representation. 
We have given the expression for the free energy of a pure gauge SU{N)/Zf^ 
model at strong coupling. Adjoint Wilson loop expectation values can be 
read out immediately for these theories. We then proceeded to show how the 
large mass expansion of theories with matter fields in the adjoint represen- 
tation yields an effective action suitable for treatment with our formula. We 
illustrate the apphcations by calculating explicitly the large N corrections 
to various quantities in the lattice version of Supersymmetric Yang-Mills. In 
all cases, to obtain a compact form, one is faced with a path resummation 
problem which at the moment we have been unable to solve. (The problem 
looks tractable, though). We also showed that our present formulas can be 
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used to prove the form of the leading large behaviour for p-gluino propa- 
gators and spectra given in Ref. [|]. We then took the case of the gaussian 
Kazakov-Migdal model to show how, similarly, one arrives at the large 
values of the masses of the p-particle bound states at strong coupling. 

From the methodological point of view, our paper might also prove of 
interest in other contexts where the large N limit is at play. Our diagram- 
matical technique is non-standard and could help to complement other well- 
known techniques as loop equations, orthogonal polynomials, master field 
equation, etc. It allows the computation of higher order corrections in l/N"^ 
keeping the Feynman diagram analogy. On the negative side, such a tech- 
nique does not look appropriate for deahng with the weak coupling (large 
5(7)) behaviour of the generating function. 

Acknowledgements 

This work has been financed by the CICYT under grant AEN97-1678. A. G- 
A wants to thank the Theory Division at CERN for its hospitality during 
the last stages of this work. The authors acknowledge useful conversations 
with E. Gabrielli and J. Vermaseren. 

Appendix A 

In this appendix we study group integrals with variables in the adjoint rep- 
resentation of the gauge group. For this purpose, as explained in Section 2, 
we will relate the matrices in the adjoint representation to their fundamental 
counterparts, and afterwards we will use known integration formulas for the 
latter. These will be first applied to the computation of single link integrals, 
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and then we will use the result to evaluate the multilink integral in Eq. 
and reduce it explicitly to the compact form in Eq. d). 

Let us first establish some notations and conventions. We write A"^ for 
the generators of SU (N) in the fundamental representation. They close the 
group algebra: 

[A", A''] = if^'X' , (27) 

where the indices a,b,c range from 1 to A^^ — 1. We choose the generators 
to be normalised as: 

Tr (A"A'') = -6"-^ . (28) 
Furthermore, the completeness relation is written as: 

Now, the group integrals we want to compute have the general form: 

j dnH{g)W^^^ . . .W-^- , (30) 

where U = Uyi{g) is the matrix corresponding to a given group element g in 
the ac^jo^JT-t representation, and dfinig) is the Haar measure on SU{N). 

The relation of U with the matrix U = Upig) corresponding to the same 
group element in the fundamental representation reads, with the conventions 
adopted above: 

= 2 Tr (f/U"t/A^) . (31) 

This allows us to rewrite (|30D as an integral for matrices in the fundamental 
representation, whose indices are contracted with a certain tensor built up 
of a product of generators. Before casting this statement into an explicit 
equality, let us introduce the following known formula for group integrals in 
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the fundamental representation 



m 



J dfiM u''^^ . . . w^^^iwy^""^ . . . (f/t)^-'"- = 

J2 S{j,a{l))6{m,a'{z))C^^{aoa') , (32) 
with n < N . 

In this expression, S„ is the permutation group of n elements, and we have 
introduced the following notation: 

n 

5(j,^(0)=nw) ' (33) 

k=l 

where la{k) is the fc-th element of cr({/i, . . . , /„}), with a G S„. 

The coefficients are class-functions on the space of permutations, i.e., 
C]^ is a map S„ ^ R fulfilling o a' o a'^) = C^ia') for any {a, a') 

(the notation a o a' stands for the composition of two permutations), so 
that the value C'^{a) only depends on the total number of cycles in the 
permutation cr. For this purpose we call Ckic) the number of fc-cycles in a, 
and C(cr) = J2k ^k{o')- With this notation, it can be shown that the following 
relation holds [|lO[] : 

6.'^= 5^C]^(a'oa)iV^('^°'^) , (34) 

creSn 

which allows to determine C^{a) recursively. Its general form turns out to 
be C]^(a) = ™, where Q{N) = N\N^ - 1) . . . {N^ - {n - if), and P{N) 
is a polynomial whose leading term is N'^^"\ Hence, at large N we get: 

C]v(c^) ^=°° -j^Scr^e + subleading terms , (35) 
where e is the identity permutation, which maps a sequence to itself. 
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If we now merge Eqs. (^Tj) and (P^), the result for (^) can be written as: 
y"d^^W''l''^..W""''" = 2" C^^ia oa)T{ay^-''-T{af^-^- , (36) 

where we have introduced the following notation: 

(37) 

an implicit summation over matrix indices being understood. These objects 
are products of traces of products of generators, and therefore invariant group 
tensors. To make this more explicit, we can introduce the one-step-forward 
cyclic permutation of m elements L^-*, defined by L^^\ai) = Oj+i for i < m 
and L^^(am) = ai, such that: 

^^l_M)ai...a^ = Tr (A"^ . . . A'^'") . (38) 

Then we notice that each cycle in a generic permutation of n elements a, as 
that appearing in (0), closes an L+ trace, so that we can write: 

CM 

T((t)"i-"" = JjT(Lf'^)''->w-'*-»(-«) , (39) 

1=1 

where ai is a sub-permutation giving the index order inside the i-th cycle, 
which is a i/j-cycle. 

Now we address the computation of Eq. (§). If we substitute for each 
single-link integration in it the result (^) we get: 

^(71 . . . 7n) = 2^ 5^ s{pr-a.b.-b. -Q [ ^ 



k=i (T(fc),5{'=)es„ 



k 
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where the index k labels the links which contribute at least one U in Eq. , 
the total number of these links being C (note that the integral over a link 
carrying no W is 1). We call the number of Ws in (^) lying on the link k; 
as we define in Section 2 the total number of Ws to be L, it is obvious that 
Ef=i = L. 

One can put the result in a more compact form by introducing a permu- 
tation ^ E S2L combining the set {a^''\a^''^; k = 1, . . . ,£}. It is important 
to notice that the set Sr of all the possible i) obtained in this way is, in 
general, only a subgroup of S2L, determined by the structure of the joint 
path 7i U . . . U 7„. Notice that the condition for Cj^f^-d) to be leading in 
is read from Eq. ( pS]) to be a'-'^^ = a'-'^-* ^ for any k; when written in terms of 
the joint permutation this condition is equivalent to i^T'dr = e, where r is 
a pairing, defined in Section 2, such that r(aj) = bi, rihi) = ai. 

Once d is fixed, we can define: 

c 

Cj,{d) ^ 2^n^^'(^^'^°^^'^) ' (41) 

k=l 

C 

'X(J)yi-C2L ^ J~j'^(^^(fc)^aia2...a„^-2'(aW)^l^2...fen;, ^ ^42) 

k=l 

after renaming the colour indices as in Section 2. It is these definitions, 
together with Eq. (|), that give the final form Eq. (|^). 



Appendix B 

In this appendix we will study the large behaviour of scalars (invariants) 
formed by contraction of invariant tensors. This type of objects appears 
in the definition of the quantities entering in Eq. (^). We will use the 
symbolic representation for tensors as given in (|37|). For notational details we 
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refer the reader to Appendix A. Computational methods for calculating these 
quantities have been devised |]14| due to their interest in different contexts. 



Our iterative technique is related to the method of Cvitanovic |]15 . 

First of all we will set up a Feynman diagram-like representation for the 
invariant tensors 7"(cr)"i - "m. For that, one uses the decomposition of the 
permutation a into cycles (see Eq. (^)). The tensor corresponding to an 
m-cycle T'^L^^)^^ - '''" is the trace of the product of m group generators, as 
shown in Eq. (^). Now, to each L^"* tensor we associate an m-legged vertex, 
where the external legs are attached to indices, oriented clockwise according 
to the orientation given by the cycle. In this way, a tensor associated to the 
permutation a is associated to C{a) vertices. Furthermore it will prove useful 
to consider vertices with no legs, associated to the trace of the identity. This 
representation is illustrated by the following: 



b„ b, 



r(L 



Tia' 



bl...bn ^ ■ /}-^^^a,(3) • ^^^r" ''''am(3) 



Tr (I) = 




The contraction of two indices of a tensor is represented in the diagrams 
by joining the corresponding legs. Given a product of invariant traces with 
contracted indices, these prescriptions allow to represent it by a diagram con- 
sisting of a certain number of vertices, joined by internal lines, and possibly 
with some external lines representing unsaturated indices of the tensor. For 



28 



the scalar quantities we are interested in, the corresponding diagram is a 
vacuum diagram with no external legs. Notice, however, that the ordering 
of the legs in a vertex is relevant, up to cyclic permutations. As explained 
in Section 2 and used later, the actual pattern of contraction of indices in 
a tensor can be labelled by a pairing p of the indices, that is a permutation 
made up of only 2- cycles. 

Our technique is based upon equations which express the result of con- 
tracting two indices in an (m + + 2)-tensor as a linear combination of 
(m + ?T,)-tensors. These equations can be obtained easily from the complete- 
ness relation of the group generators Eq. (pQ]). We will refer to the basic 
equations as fusion and fission rules, and they are given respectively by: 

q-^y^im+l)ybi...b„iq-^y^{n+l)yci...c„ _ ^ q- ^^(m+n)^bi...b,nCi...c„ _ 

a 

q-^^^{m+n+2)ybi...b„,aci...c, 



-T(L 

2 ^ 



{m)^^6l...^)„^n {n)-vci...c„ 



{m+n)^bi...b^ci...Cr, 



2N 



t(l!^ ) 



(44) 



In our diagrammatic representation the fusion and fission rules can be 
depicted as: 



It is easy to see that the repeated apphcation of Eqs. P3| , |iB| ) allows the 
evaluation of all scalars formed by contracting the indices of an arbitrary 
tensor. Our goal will be to extract information on the dependence from 
these equations. As an example and for later use, we give in Table 1 the 
value and associated diagram for all independent tensor contractions up to 
six indices. The corresponding diagrams have K = 1 connected components 
and up to L = 3 internal lines. 

If the number of indices of a tensor is odd, making all possible contractions 
would leave an unsaturated index. Applying the fusion-fission rules would 
reduce it to a linear combination of tensors with one index. But the only 
one-index tensor Tr (A") vanishes, and hence the contribution of a diagram 
with a single external leg vanishes too. This allows to prove a result which 
will be needed in what follows: the scalar corresponding to a one-particle 
reducible diagram vanishes; by applying fusion and fission rules iteratively 
in each of those two sides of the diagram connected by the reducible line we 
would reach the form B{N) Tr (A'^)Tr (A"), which is zero q.e.d. 

Now we will prove a result concerning the A^ dependence of the scalars 
formed by contraction of a given tensor, expressed in terms of the topological 
properties of the associated diagram. The result is that the leading A^ be- 
haviour of such a scalar is given by m2~^N'^^~^~^^ , with m being either or 
1. It is easy to check this property directly for tensors with few indices. The 
tensor contractions given in Table 1 and their multiplication allow to check 
it up to L < 3. The general proof will be done by induction, making use of 
the fusion and fission rules. Hence, we will start by assuming that the result 
has been proven for all scalars obtained after the contraction of L pairs of 
indices, and end up proving that the result is also true for L + 1 contractions. 

Consider then a scalar formed by performing L + 1 contractions of the 
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indices of a 2L + 2 invariant tensor and its associated one-particle irreducible 
diagram. Let say it has K connected components, V vertices and L + 1 
internal lines. Applying the fusion rule its contribution is given as a linear 
combination of the contribution of two diagrams with L internal lines, with 
coefficients 1/2 and —1/2N. The first diagram has L internal lines, V — 1 
vertices and K connected components; thus, its leading term is of the form 
mi2~^A^^^~^+^+^. The second has L internal lines, V vertices and K con- 
nected components^, so that the leading contribution is m22~^N'^^^^^^^^ , 
which is subleading with respect to the previous one. In summary, the be- 
haviour of the initial scalar is given by rn2~^^~^^^ N'^^~'^'^^^'^^\ as predicted 
by our result, with m = mi. 

Now let us turn to the fission rule. As the original diagram has K con- 
nected components, V vertices and L + 1 internal lines, then the two sub- 
diagrams in the decomposition have {K', V + 1,L) and {K, V, L) topological 
numbers respectively, with K' being either K oi K + 1. The leading term 
of the original diagram is thus of the required form m2~^^^'^'> N'^^~'^^'^^^^\ 
with m = if K' = K and m = mi if i^' = + 1. This completes the proof. 

The previous result is used in Section 2, to prove that the coefficients 
Dj^l'd, p), which are associated to a given diagram, have the mentioned large 

behaviour. This, combined with the leading behaviour of the coefficients 
Cm, leads to the conclusion that the large A^ dominant behaviour of the 
integrals S (Eq. (^)) can only come from diagrams made out of connected 
components with two vertices. From our previous results we also know that 
the prefactor of the leading term is an integer counting the number of dia- 
grams giving a non-vanishing maximal leading behaviour. For that purpose 

■'it cannot have K + 1 since this would contradict the hypothesis that the original 
diagram is one-particle irreducible. 
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we need to complement our previous result by characterising those diagrams 
which have a maximal large N behaviour {m — 1). Fortunately, we only need 
to do this for diagrams with two vertices V = 2 (and no self-contractions), 
making up one of the connected components of a possible maximal diagram. 

In the remaining of this Appendix, we will prove the statement made 
in Section 2, namely that the two-vertex connected (and with no self-con- 
tractions) diagrams having maximal leading behaviour (N^) are the ones 
satisfying the identity 't}frdp — e, where is the permutation giving the 
vertex structure and p the permutation specifying the contractions. By the 
form of the diagram we know that -& is made of 2 cycles of length L, and p 
of L cycles of length 2. Let us see what this means in diagrammatic terms. 
Take the legs of one vertex and label one as ai, then label the subsequent legs 
Oj+i = '&{ai). This means that the legs are ordered in a clockwise fashion. 
Now call bi = p{a.i) the indices of the other vertex (remember there are no self- 
contractions). The condition dpdp = e is equivalent to having d(hi) = bi-i. 
This means that in the other vertex the ordering of the indices is counter- 
clockwise. The corresponding diagram can then be drawn by joining the legs 
with lines which do not cross: the diagram is planar. 

Using this graphical characterisation, it is easy to show that a diagram 
coming from e can never have a leading term, of the form studied 

above, with m = 1. On the other hand, an exact formula can be derived for 
planar diagrams by iterating the fusion rule; the result with L > 2 lines is: 

(^)'O-(r-^) . (*^) 

and obviously gives leading behaviour with m = 1. 

To prove that any nonplanar diagram is necessarily not leading, let us 
depict a generic nonplanar diagram with L — s + 2 lines in the following 
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form (which is always attainable), 




where the indexed legs are completely contracted between themselves 
(remember that no self-contractions are allowed), and s — m + n— p + q> 1. 
The nonplanarity condition enters as a constraint: both a — c and a — d 
contractions exist. 

If a fusion is applied, and then a fission to the resulting first term in the 
rhs, one gets: 




Now, the nonplanarity constraint prevents the first term in the rhs to have 
more than one connected component. It is then immediate, after applying 
the rule for the leading term to each diagram in the rhs, that the leading 
behaviour for the nonplanar diagram goes at most like N^"^. 
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Table 1: Explicit results for all topologically inequivalent, nonvanishing, com- 
pletely connected diagrams corresponding to contractions of invariant traces 
with up to 3 internal lines are shown. Disconnected diagrams with L < 3 
can be built as products of the ones shown here. 









= Tr(A''A") 

= Tr (A"A''A''A") 

= Tr (A"A''A"A'') 

= Tr (A"A^)Tr (A'^A^) 

= Tr (A"A''A'=A"A''A'^) 

= Tr (A"A"A''A''A^A'^) 

= Tr (A^A^A^A^A^A'^) 

= Tr (A^A^'A^A'^A^'A'^) 



2 

In 

AN 

N^-1 
4 

{N^ - 1)3 

(A^^ - If 

(iV2 - 1)2 



N^-1 
8m 



.a^^>^c^^.Mc^.^a^ _ (A^^ - l)(iV^ - 2) 



= Tr (A"A''A")Tr (A"A'A° 



8N 
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AN 



Tr (A^A-^A^A^jTr (A''A'=) = 



(iV2 - 1) 

8N~ 



Tr (A"A''A"A")Tr (A^'A") = 



SN 



y (A"A'')Tr (A''A")Tr (A"A") 
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